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ABSTRACT 

The  general  equations  of  surface  ship  motion  are 
developed  and  standardized  for  simulation  in  digital 
computer.   Digital  simulations  of  the  dynamics  of  the 
surface  ship  in  three  degrees  of  freedom  are  done 
with  and  without  non-linear  and  cross-coupling  terms. 


TABLE  OF  CONTENTS 

I.  EQUATIONS  OF  SURFACE  SHIP  MOTION  6 

II.  DIGITAL  COMPUTER  SIMULATION  15 

III.  STUDY  OF  SHIP  "D"  PERFORMANCE 27 

A.  NON-LINEAR  TERMS  NOT  INCLUDED  30 

B.  INCLUDED  NON-LINEAR  TERMS  44 

IV.  CONCLUSIONS 63 

COMPUTER  PROGRAM  32 

1.  SURFACE  SHIP  3  DEGREES  OF  FREEDOM 

WITHOUT  NON-LINEAR  TERMS  32 

2.  SURFACE  SHIP  3  DEGREES  OF  FREEDOM 

WITH  NON-LINEAR  TERMS  46 

3.  SUBROUTINE  TO  CALL  DRAW 61 

LIST  OF  REFERENCES 64 

INITIAL  DISTRIBUTION  LIST  65 

FORM  DD  1473 66 


LIST  OF  TABLES 


I.  Hydrodynamic  coefficients  and  constants 

of  ship  "D"  for  linear  terms 30 

II.  Hydrodynamic  coefficients  of  ship  "D" 

for  non-linear  terms  45 


ACKNOWLEDGEMENT 

The  author  wishes  to  express  his  sincere  appreciation 
to  Professor  MILTON  L.  WILCOX  for  the  guidance,  assist- 
ance and  continuous  encouragement  which  he  provided  during 
the  pursuit  of  this  study.   The  author  would  also  like  to 
express  his  appreciation  to  Dr.  GEORGE  J.  THALER  for  his 
valuable  assistance  in  furnishing  the  data  for  simulation. 


I,   EQUATIONS  OF  SURFACE  SHIP  MOTION 

A  moving  ship  is  a  body  with  six  degrees  of  freedom. 
These  degrees  of  freedom  are  generally  chosen  as  follows: 

a.  Linear  displacements  along  the  three  axes 
through  the  center  of  gravity. 

a.l  Surge  —  along  X  axes 
a. 2  Sway  —  along  Y  axes 
a. 3   Heave  —  along  Z  axes 

b.  Rotations  around  the  three  axes  through  the 
center  of  gravity. 

b.l   Roll   —  around  X  axes 

b.2   Pitch  —  around  Y  axes 

b.3   Yaw    --  around  Z  axes 
Further  reduction  in  the  complex  nature  of  the  equa- 
tions can  be  brought  about  by  choosing  an  orthogonal 
axis  system  parallel  to  the  principal  axes  of  inertia  so 
as  to  eliminate  products  of  inertia  in  the  motion  equa- 
tions.  For  practically  all  ocean  vehicles,  with  extreme- 
ly few  exceptions,  a  longitudinal  axis  (X  axis)  in  the 
centerline  plane,  a  downward  (toward  keel)  axis  (Z  axis) 
perpendicular  to  the  X  axis  in  the  centerline  plane, 
and  a  transverse  axis  (Y  axis)  perpendicular  to  the 
centerline  plane  satisfies  this  requirement. 
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Figure  1.   Surface  ship  in  six  degrees  of  freedom. 


For  the  exceptional  vehicle  which  has  a  very  peculiar 
and  significantly  large  asymmetrical  mass  distribution, 
it  is  necessary  to  include  the  products  of  inertia. 

The  X,  Y,  Z  axes  form  an  orthogonal  right  hand  system 
of  axes  fixed  in  the  vehicle.   The  axes  and  the  associated 
components  of  the  pertinent  physical  quantities  are 
defined  below: 

The  longitudinal  X  axis  (in  the  plane  of  symmetry)  is 
positive  in  the  forward  direction,  usually  parallel  to 
the  keel  or  calm  water  line.   If  the  upper  and  lower  halves 
of  the  body  are  symmetrical,  then  the  axis  is  the  inter- 
section of  the  two  planes  of  symmetry. 

The  Y  axis  is  the  transverse  axis  perpendicular  to  the 
plane  of  symmetry  and  positive  to  the  starboard. 

The  Z  axis  or  downward  axis  in  the  plane  of  symmetry 
(X,Z)  is  perpendicular  to  the  X  axis  and  positive  downward 
towards  the  keel. 


>\    S\  >N 


i,j,k  unit  vectors  along  the  X,Y,Z  axis  respectively. 
R  x,y,z  vector  distance  of  a  point  from  the  origin  0, 
and  the  corresponding  components  along  the  X,Y  and  Z  axes. 

R  =  ixG  +  jyG  +  kzG 

U  u,v,w  velocity  or  the  origin  0  (on  the  body)  and 
the  corresponding  components  along  the  X,Y  and  Z  axis. 

/\  /\        y\ 

U=iu+jv+  kw 
»  P/q/r  angular  velocity  of  the  body  about  the  origin 
and  the  corresponding  components  about  the  axes. 

v       s\  s\  y\ 

u   =  ip  +  jn  +  kr 

The  moments  of  inertia  of  the  body  about  the  X,Y  and 

Z  axes  respectively  1,1,1. 
r       J      x   y   z 

F  X,Y,Z,  force  acting  on  the  body  and  the  correspond- 
ing components  along  the  axes. 

y\  /\         s\ 

F  =  iX  +  jY  +  kZ 
— > 
m  K,M,N  Moments  acting  about  the  axes. 

m  =  lK  +  jM  +  kN 

Newton's  law  of  motion  for  a  rigid  body  can  be  written 

as  two  equations ,  one  a  force  equation  and  the  second  a 

moment  equation  provided  an  origin  is  taken  at  the  center 

of  gravity  and  the  axis  system  is  fixed  in  space.   The 

equations  are 


F  =  -T77  (momentum)  •  prp  (m  U_) 


m  =  gx-  (angular  momentum)   =  -=-r-(Ifi) 


where  the  subscript  G  refers  to  an  origin  at  the  center 
of  gravity  and  m  is  the  mass  of  the  body.   For  a  mass 
essentially  constant  in  time 

For  an  origin  not  at  the  center  of  gravity  of  the 

body  and  in  a  system  of  axes  fixed  in  and  moving  with 

the  vehicle. 

U_  =  U   +  S  x  R_ 
G     a         G 

where  U   is  the  velocity  of  the  origin  in  space.   However, 
a 

since  the  origin  is  on  the  surface  of  the  earth  and  the 

earth  rotates ,  then 

->    ->-->-    -> 
U   =  U  +  ft   x  R, 
a         e     b 

where  U  is  the  geographical  velocity  of  the  body,  ft   is 

the  angular  velocity  of  the  earth,  and  R,  is  the  radius 

vector  from  earth's  center  to  the  vehicle.   The  force 

equation  becomes: 

F    =   m   -^(U    +    ft      x    R,     +    ft    x    R    ) 


±  d     ,->        -*■        -*■    .        ,4-        -*■-*■  4-    , 

F   =   m^r-(U   +   ft    x   R)+m[ft    x   R,+ft    x  R,] 

=   m   orr(U    +   ft'   x   R_)+m[&    x   U   +fi   x  ft*   x   R,  ] 
dt                         G             e                e        e        d 

since   ft    =    0    and   R,     =   U   +   ft      x   R, 

the    term  mft    xU      is    the    coriolis    force    and   mfi*      x   ft  x    R, 
e                                                                               e           e  b 

is    the    centripetal   acceleration   due    to   rotation   of  the 


earth.   These  two  terms  are  negligibly  small  when  com- 
pared with  the  other  forces,  then 


d_ 

dt 


—  (U  +  ft  x  R_) 


Finding  the  derivatives  of  unit  vectors  (change  in 

direction) 

where     di  =  -kd  0      di  =  jdt       di  =  0 

dj  =  0         dj  =  -id^     dj  =  kd  |) 

dk  =  id  9       dk  =  0        dk  =  -jdj> 

Adding  the  contributions 

di        ?    n        ?dijj         ~d9 
dt  Jdt  dt 

«i  =  _£a  +  j-o  +  k|| 

dt        J  dt 


dt 

dk 
dt 


?  d_e_ 

dt 


=  i  tt  -  jarJ-  +  k-0 


^d4> 
'dt 


.  /\  /v  /\ 

u    =    ip    +    jq    +   kr 


dcf> 


dt    ' 


q  = 


de 

dt    ' 


r    =    fflL 
dt 


di 

dt 


=   ft 


x    1    = 


P 
1 


q 


k 

r 


0         0 


du 

dt 


U  =   iu   +  jv   +  kw 

*G   =    *XG+  ^YG+    *ZG 

?  ,  .Cdi        "•  di 

=   in  +   u(j-r-  +    jv  +  v=f-  +   kw   + 


dt 


dt 


w 


dk 
dt 


=   iu  +  u(5||-^)    +jv+v,-2^+^, 


■dt 


■dt 


+        kW       +       W(l^-        Dg£) 


=    iu+    u(jr-   kq)+    jv+   v  (-ir+   kp)+kw   +   w(iq    -    jp) 
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3 


x\   = 


•s.  /N  /S 

i         j         k 

p      q      r 


XG   YG    ZG 


i(qZr-    rYr)-    j  (pZ    -    rX_)+    k  (pY    -qX,J 


—  ($^xRG)    =    i(ZQq-YGr)+    j  (ZGqr-YGr2 )  +   k(YQqr-ZGq2) 

S\  /N  #-j  /S.  *-v 

+  j  (XGr  -  ZGp)+  k(XGrp  -  ZGp  )+  i(ZGrp  -  XQr  ) 

+    k(YGp    -    XGq)+    i(YGpq    -    XQq    ) +    j  (XQpq    -    Y^) 
from        F  =   ix   +   jY   +  kZ      and     F  =  m[jS-  +   ^kr(^   x   RQ) 

Rewriting  and  grouping  all  i  terms  equal  to  X,  all  j 
terms  equal  to  Y  and  all  k  terms  equal  to  Z  yield 

X  =  m(u-vr+wg+ZGq  +  rpZG~r2XG+  pqYQ  -  q2XQ) 

2   2  •  • 

=  m(u+wg-vr-XG (r  +q  )  +  YG(pq-r)  +  ZG(rp+q) 


Y  =  m(v+ur-wp+XG(pq+r)  -  YQ(p2+r2)  +  ZQ(pq-r) 
Z  =  m(w+vp-uq+XG(rp-q)  +  YG(qr+p)  -  ZQ  (p2+  q2) 


(1) 


From 


G  =  -^(angular  momentum) 


m„  = 


=  5t(i:cx  P  +  ^y  q  +  kIz  r) 


G  indicates  an  origin  at  the  center  of  gravity 

JX  "  *X  -  »  (YG2  +  ZG2) 
G 

iYq=  iy  -  »  (zG2  +  xG2) 

h    -  h    ~   m    (XG2  +  YG2) 

M  =  MG  +  RQ  x  F    or   MQ  =  M  -  R   x  F 
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After  manipulating  and  using  the  results  for  the 
derivatives  of  unit  vectors  (same  as  above)  expressions 
for  K,  H  and  N  are  obtained. 
K  =  Ixp  +(IZ-  Iy)qr+m[YG(w+  pv  -  qu) -  ZQ(v+ru-pw)] 

M  =  Iyq  +(IX-  Iz)rp+mtZG(u+  qw  -  rv) -  XQ (w+pv-qu) ]    (2) 

N  =  Izr  +(IY-  Ix)pq+m[XG(v+  ru  -  pw) -  YG(u+qw-rv)] 

The  terms  (qw-rv) ,  (ru-pw)  and  (pv-qu)  are  gyroscopic 
effects . 

The  relationship  for  forces  and  moments  can  be  ex- 
pressed 

Force) 

)  =  f (properties  of  body,  properties  of 
Moment)      motion,  properties  of  fluid) 

For  a  particular  ship,  in  a  given  fluid  with  no 
excitation  force  -  so 

*  ) 

_^  )  =  f  (properties  of  motion) 
M  ) 


•    m  • 


=  f  (Xo,Yo#Zo,(j),e  ,^,u,v,w,p,q,r,u,v,w, 
p,q,r,6,6\S") 
The  Taylor  series  which  has  the  following  form  may 
now  be  applied  to  linearize  the  equations  about  an 
operating  point  X 

f(X)  =  f  (X  )  +  (X-X  )%-   f(X  )+(X-X  )2  ^-,  f(X  )  + 

dX 

Apply  this  to    f(X,Y,Z) 
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For  the  case  of  f (properties  of  motion)  let  (X-X  ) 

=  AX  ,  (Y-Y  )  =  AY   and  (Z-Z  )  =  AZ 
o      o      o         o      o 

terms  second  order  and  higher  are  neglected  for  small 

perturbations . 

From  X  equation  the  linear  terms  are  obtained: 

X   =   f(— )0+   AX0(||  )+    4Y0<||  )+    4Z0(«    )  +    ...AV(||)+    ... 

o         o         o 

The  defining  relations  are: 

Wo  =  Wu=u         =  Xu 
o 


(fe  =  X 

cw  w=w  w 

o 

Aw  =  (w-w  )=w,    w   =  0 
o        o 

Au  =  (u-u  ) 

The  force  equations  then  become 

X    =    X      +    Xv    X      +    X.,    Y      + XQ0    +     X    A 

0X0Y0  0  uu 

o  o 

Y  =   Y      +YVX      +Y^Y      + YQ6    +     Y    A  (3) 

oXoYo  6  uu 

o  o 

Z=Z      +    Zv    X      +ZVY      +....Z„6    +     ....Z    A 
oXoYo  9  uu 

o  o 

A  similar  derivation  can  be  done  for  the  K,M  and  N  equa- 
tions.  The  preceding  X,Y  and  Z  equations  may  now  be 
equated  to  the  linearized  X,Y  and  Z  (equations  (1)),  e.g., 
for  the  Y  equation  without  roll,  pitch,  and  the  center 
of  gravity  at  X  =  0,  Y  =  0  and  Z  =  0  gives  the  linearized 
equation . 

Y  =  m(v  +  ur) 
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then 

Y„  X  +  Y„  Y  +  Y.iH-  Y  u+  Y  v+  Y   r+Y-u+   Y«v+   Y-r  =  m(v+ru) 
XoYoiliu  v  ru  v  r 

o  o  T 

Expressions  for  X,Z  and  K,M  and  N  can  be  determined  in 

a  similar  procedure. 

In  order  to  obtain  the  equations  in  a  non-dimensional 

form  some  definitions  will  be  given,  and  applied  to  the  Y 

force  equations  as  an  example  of  the  process. 

Froude  number  =  


/gT 


i        i        i         u,v,w 
/g/T 


f         t(/g/£) 
X,Y 

pg£- 


X',Y' ,Z'     =    X'Y^Z 


^X    AY    '^Z  1      TT2„2 

y   pU    £ 

After  replacing  and  adding  the  effect  of  waves,  the 
Y  equation  becomes: 

v'  +  r'u'  =  \   U2, (Cv   v'  +  Cv   p'  +  Cv   r '+  r        fir) 

v'       Yp'       Yr'      Y6r 

+  (Y^.p1  +  Y^  v»  +  Y-,r'  +  Y  ») 
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II.   DIGITAL  COMPUTER  SIMULATION 

The  six  equations  of  motion  after  rearranging  by  plac- 
ing the  second  order  terms  to  the  left  and  the  rest  on  the 
right  side  become: 

aaA+baB+caC+dab+eaE+faF  =  -  (a.ja,  A+a2a2A+b,  a.jB+b2a2B 

+  c,a,C+c2a  C+d-ja-jD+d^^ 
+  e1a1E+e2a2E+f1a1F+f2a2F) 
+  IF1 
abA+bbB+cbC+dbD+ebE+fbF  =  -  {a^^A+a^^+b^b-jB+b^^B 

+  c1b,C+c  b  C+d,b1D+d2b2D 
+  e1b1E+e2b2E+f1b1F+f2b2F) 
+  IF2 
acA+bcB+ccC+dcD+ecE+fcF  =  -  (a,  c,  A+a^c-A+b..  c.  B+b^c^B 

+  c1c,C+c  c  C+d1c  D+d  c  D 
+  e1c1E+e2c2E+f1c1F+f2c2F) 
+  IF3 
adA+bdB+cdC+ddD+e'dE+f  dF  =  -  (a,d  A+a2d  A+b,d  B+b  d  B 

+  c1d1C+c2d2C+d1d  D+d  d  D 
+  e1d1E+e2d2E+f1d1F+f  d  F) 
+  IF4 
aeA+beB+ceC+deD+eeE+feF  =  -  (a^-jA+a  e  A+b  e,B+b  e  B 

+  c,e1C+c2e2C+d1e1D+d2e  D 
+  e1e1E+e2e2E+f1e,F+f2e2F) 
+  IF5 
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afA+bfB+cfC+dfD+efE+ffF  =  - (a,f 1A+a2f 2A+hif iB+b2f 2B 

+  e1f1E+e2f2E+f1f1P+f2f2P) 

+  IF6 
where  A=u,  A=u,  B=v,  B=v,  C=W,  C=W,  D=p,  D=p 

E=qf  E=q,  F=r ,  F=r ,  terms  IF  include  all  non-linear 

terms  such  as  wave  force,  wind,  rudder  deflection etc. 

In  the  six  equations,  the  right  can  be  set  equal  to 


II'  X2' 


■Ifi  respectively,  thus: 


I,  =  - (a, a, A+a„a?A+b, a, B+ 
I2  =  -  (a-jb,  A+a^b-A+b,  b-,B+ 


I-  =  -  (a  c  A+a  c  A+b,c1B+ 

I.  =  -  (a,d,A+a„d^A+b.d,B+ 
4        11    2    /.  j.  j. 

Ic  =  -  (a, e, A+a-e-A+b, e, B+ 

Ic  =    -  (a, f ,A+a0f„A+b. f , B+ 

D  11      2      2  11 


f2a2F 
f2b2F 


f2C2F 
f2d2P 

f2e2F 
f2f2F 


+  IF1 
+  IF2 

+  IF3 
+  IF4 
+  IF5 
+  IF6 


the  equations  then  have  the  form  that  follows , 

aaA  +  baB  +  caC  +  daD  +  eaE  +  faF  =  I 

abA  +  bbB  +  cbC  +'  dbD  +  ebE  +  fbF  =  l? 

acA  +  bcB  +  ccC  +  dcD  +  ecE  +  fcF  =  I 

adA  +  bdB  +  cdC  +  ddD  +  edE  +  fdF  =  I. 

4 

aeA  +  beB  +  ceC  +  deD  +  eeE  +  feF  =  I 
afA  +  bfB  +  cfC  +  dfD  +  efE  +  ffF  =  I, 
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expressing  in  matrix  form: 


aa 

ba 

ca 

da 

ea 

fa 

ab 

bb 

cb 

db 

eb 

fb 

ac 

be 

cc 

dc 

ec 

fc 

ad 

bd 

cd 

dd 

ed 

fd 

ae 

be 

ce 

de 

ee 

fe 

af 

bf 

cf 

df 

ef 

ff 

.. 

p 

A 

xi 

B 

T2 

C 

X3 

D 

J4 

E 

J5 

F 

h 

Apply  Cramer's  rule  to  solve  for  A,  B 


F  in 


terms  of  I,  --  I,- 


A  = 


h 

ba 

ca 

da 

ea 

fa 

l2 

bb 

cb 

db 

eb 

fb 

J3 

be 

cc 

dc 

ec 

fc 

Z4 

bd 

cd 

dd 

ed 

fd 

X5 

be 

ce 

de 

ee 

fe 

h 

bf 

cf 

df 

ef 

ff 

aa 

ba 

ca 

da 

ea 

fa 

ab 

bb 

cb 

db 

eb 

fb 

ac 

be  • 

cc 

dc 

ec 

fc 

ad 

bd 

cd 

dd 

ed 

fd 

ae 

be 

ce 

de 

ee 

fe 

af 

bf 

cf 

df 

ef 

ff 

define  the  denominator  determinant  _  A  and  for  the 
nominator  let  cof actor  of  I,  =  cof.  aa 


cof actor  of  I_  _  cof.  ab, 


=  cof.  af  equations  becomes: 


•and  coefactor  of  I. 


A  = 


(cof  .  aal,  +cof  .abl^+cof  .  ad_,+cof  .  adl . +cof  .aeIr+cof  .  af  I,.) 
1         2         3         4         5         6 
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In  the  same  way  solve  for  B,  C,  

(cof  bal.  cof  bbl-  cof  bd3  cof  bdl.  cof  bel.  cof  bflg) 


B  = 


A 


(cof  cal,  cof  cbl~  cof  ccl-,  cof  cdl.  cof  eel,-  cof  cfl_) 
" 1 2 3 4 5 6 

L  ~  A 

(cof  dal,  cof  dbl~  cof  dcl0  cof  ddl.  cof  del-  cof  dfl..) 
U  ;  A 

(cof  eal,  cof  ebl„  cof  eel-,  cof  edl.  cof  eelc  cof  efl,) 
„           i        2                    o                    4         o         b 
E  _   _ _ 

(cof  fal,  cof  fbl~  cof  fcl_.  cof  fdl.  cof  fel-  cof  ffl^) 
t-,  1  2  J  4  b  b 

F  =   ! 


Then  the  value  of  A,  A,  B,  B F,  F   by  integration 

such  that 

A  =  5Sl  =  f  £h  a    -  \  *-A 

A    dt    J   dt2  '   A  ~  J  dt 


A  block  diagram  to  compute  all  of  the  variables  in  the 
set  of  equations  is  presented  in  Fig.  2. 

In  the  computer  program  that  is  used  for  simulation 
all  six  equations  for  six  degrees  of  freedom  are  used,  but 
are  interested  in  less  than  six  degrees  of  freedom.   The 
same  program  can  be  used  by  setting  the  coupling  terms  of 
non-used  equations  equal  to  zero  and  one  in  terms  of  princi- 
pal diagonal,  e.g.  only  three  degrees  of  freedom  are  used 
in  this  study,  surge,  sway  and  yaw,  then  all  coupling  terms 
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are  set  equal  to  zero  and  unused  terms  on  the  principal 
diagonal  equal  to  one,  for  example  : 


aa 

ba 

0 

0 

0 

ab 

bb 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

af 

bf 

0 

0 

0 

fa 

•  • 

A 

Ji 

fb 

•  • 

B 

X2 

0 
0 

•  • 

C 

•  « 

D 

= 

X3 
J4 

0 

•  • 

E 

h 

ff 

•  • 

F 

Z6 

and  the  left  side  of  the  unused  equations  are  set  equal 
to  zero. 

With  this  program  non-linear  terms  can  be  added  such 
as,  rudder  deflection  of  waves  and  wind,  etc.,  which  will 
be  done  by  adding  all  of  these  whose  sum  is  IF   e.g. 
IF,  =  KA,x  Dr+KA2x  Ds+KA3x  Db+NA 

where  Dr,  Ds  and  Db  are  rudder  deflection,  canard  de- 
flection   etc.   NA  is  the  sum  of  all  non-linear 

terms  that  effect  the  surge  equation  (X  equation) . 

The  program  that  will  be  used  for  solving  these  equa- 
tions is  the  "Continuous  Systems  Modeling  Program"  (CSMP) 
[Ref.  3]  in  which  all  constants  are  declared  in  the  first 
section  and  then  set  the  value  of  matrix  for  aa,  ab,  ac 

and  so  on  (in  program  AAA  is  used  for  aa,  AAB  for 

ab  AFF  for  f f ) .   In  the  initial  section  values  of 

the  COFACTORS  aa,  ba....are  determined.   All  of  the  CO- 
FACTORS  and  the  subprogram  VALUE  is  used  to  compute.  This 
subprogram  finds  the  value  of  the  determinant  of  the 
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matrix.   For  all  of  the  COFACTOR  terms  the  element  is  set 
equal  to  one  and  the  rest  of  the  elements  in  that  row  and 
column  are  set  equal  to  zero.   For  example,  to  find  the 
value  of  COF.aa  the  following  array  is  obtained: 

10    0    0    0    0 

0    bb   cb   db   eb   fb 
COF  aa  =      0    be   cc   dc   ec   fc 

0    bd   cd   dd   ed   fd 

0    be   ce   de   ee   fe 

0    bf   cf   df   ef   ff 
(In  the  computer  program  BAA  is  used  for  a,a. ,  GAA  for 
a~a_...).   After  the  value  of  A  and  all  cof actors  are 
determined,  the  dynamic  section  is  used  to  determine  BAA, 

BAB  GAA,  GAB  (if  those  terms  contain  variables) 

In  the  dynamic  section  all  variables  that  are  functions 
of  time  are  determined.   The  defining  relations  of  the 
variables  are  also  included  in  the  dynamics  section,  i.e. 
UDOT  =  ADDOT  (U=A)  ,  U=ADOT etc.  XH,YH,ZH  are  deter- 
mined and  are  the  vector  terms,  X,Y,Z  whose  origin  is  fixed 
on  the  earth  (relative  to  the  earth) . 
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III.   STUDY  OF  SHIP  "D"  PERFORMANCE 

In  this  section  the  computer  will  be  used  to  solve 
the  equation  of  motion  describing  ship  "D".   The  hydro- 
dynamic  coefficients  and  constants  that  were  obtained 
from  NSRDC  (NAVAL  SHIP  RESEARCH  AND  DEVELOPMENT  CENTER) 
[Ref.  6]  for  this  study  will  concern  only  three  degrees 
of  freedom  such  that  SURGE,  SWAY,  and  YAW. 

Equations  of  motion  of  the  ship  (nondimensional) 
Axial  Force 


m(u-vr+wg)  =  §-  %  (X  g  +  X   r  +  X   rp) 
^2     aq^  rr       rp  tft 


gg 


+  %■  (X-u  +  X   vr  +  X   wg) 
2    u     vr      wg  ^ 

+  ^-(X   v2) 
2%      vv 

+  %^   u2(Xx  .  6r2+Xx  .  6s2+Xx,  x,  6b2) 
2 Z  6r<Sr      6s6s      Sb6b 

+  2T  (Alu2  +  A2U*Uc  +  A3Uc2)  - 


Lateral  Force 


m(v+ur-wp)  =  §-  £  (y.p  +  Y«r  +  Y   pq) 
c  2  p     r     pq 

+  £-  (Y   wp+Y  |  iv|r|+uY  r+Y-v+uY  p) 

2   wp  r  v|r|  '  '    r   v     p^ 

+  |  (UYVV  +  Ywvwv  +  Y|v|v|v|v) 

Yawing  Moment 
Izr  +  (Iy-Ix)pq  =  |  (N^r  +  N^p  +  NpgPq) 

+  ^(N-v+uN  p+uN  r+N  wp+N,  ,  Ivlr) 
21     v  p^   r   wp  -   Ivlr11 
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+  £—_-(uN  v+N   wv+N,   ,   Ivlv) 
2„ 2    v    wv     Ivlv11 

+  £-»■  u2(N.  fir) 
2£2      6r 

p,  the  density  of  fluid  is  taken  as  2  and  the  terms 
including  w,  p.  q  (heave ,  roll  and  pitch)  are  set  equal 
to  zero. 

The  nonlinear  terms  such  as  the  squared  terms  and 
product  terms  of  v  and  r  are  omitted  initially.   This  is 
in  agreement  with  the  small  perturbation  theory. 

After  rearranging,  the  equations  become 

2 


(X 


u 


\  *  „  ,    2    u 

m)  u  =   -X,    .    or      t; — 
firfir  I 


(Y-    -  m)v  +  Y'r  =   -  t  V  v   -  uY  r   -  ?-  Y.    fir 
v  r  i     v  r  I        fir 

vr  +    (N.-Iz)r   =    -   -2Nvv-   -2Nrr   -   -y  N^fir 

Set  the  left  side  of  the  equations  equal  to  I.,  , Ifi, 

then  the  matrix  equation  becomes 


X-  -m) 
u 

0 

0- 

0 

0 

(Y--m) 

V 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

N-/1 

0 

0 

0 

0 

A 

H 

0 

Y 
r 

B 

J2 

0 

0 

C 

h 

0 

0 

D 

X4 

1 

0 

E 

X5 

0 

(N--I    ) 
r      z 

F 

h 

where  u  =  A,v  =  B/  r  =  F  . 
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Also  set  the  right  side  of  the  equations  equal  to  I, 
I 


0 


0 


u. 


o  H. N    o 

%    v 


0    0 


0     0 


0 


0    ^Y     0     0     0    uY 
%   V  r 


0     0     0     0     0 


u„ 


0     0    £N 


A 
B 
C 
D 
E 

* 

F 


+ 


All  Zero 


where   u  =  A,  v=B,  r=F,  and  IFl  =  -X 


A 

IFl 

B 

IF2 

C 
D 

+ 

IF3 

IF4 

E 

IF5 

F 

IF6 

6r 


2  u* 


IF2 


6r<5r 

2 

u 


-Y.  6r 

or    1 


IF3 
IF4 
IF5 
IF6 


u 


-N.  6r 
6r    %2 
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A.   NONLINEAR  TERMS  NOT  INCLUDED 

TABLE  I 

Hydrodynamic  Coefficients  and  Constants  of  Ship  "D" 
for  Linear  Terms  (Non-dimensional) 

m  =   0.0045 

Iz  =   0.0003 

Nr  =   0.0012 

Nr  =  -0.0002 

Nv  =  -0.0012 

Nv  =  -0.0001 

Xu  =  -0.00036 

Yv  =  -0.0025 

Yv  =  -0.0063 

Yr  =   0.004 

Xx  x   =  -0.0011 
6ror 

Yx     =   0.0019 

6r 

Nx     =  -0.00084 
or 

6s  and  6b  equal  zero 

In  the  computer  program,  set  all  coefficients  in 

section  1  and  set  AAA  =  Xu-m  ,  AAB  =  0 AFF  =  Nr-Iz. 

After  this,  use  subprogram  value  to  find  the  determinant 

and  coefficient  of  AA,  BB  and  then  set  BBB  =  r-  Yv  , 

BFB  =  uYr BFF  =  J  Nr  . 

*  2 

IF1  =  KA16r   where  KA1  =  -Xx  -  6r  £— 

6r6r    I 

IF2  =  KB16r   where  KB2  =  -Yx  x? / 1 

6r    ' 

IF3  =  KF16r   where  KFl  =  -Nx   u2/£2 

6r 
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Following  is  the  program  that  used  CSMP  to  determine 
the  turning  circles  for  rudder  angles  of  15  (-0.2619  rad.), 
25°(-0.4365  rad.)  and  35°(-0.6111  rad.).   Result  of  this 
study  is  presented  in  Fig.  3.   The  turning  rate  as  a  func- 
tion of  time  is  interested  and  the  results  of  this  analysis 
are  presented  in  Fig.  4.   The  ships  turning  performance 
expressed  in  transfer  ship  lengths  as  time  is  shown  in 
Fig.  5  and  the  heading  angle  as  a  function  of  time  is 
given  in  Fig.  6.   Fig.  7  shows  the  results  of  the  zig-zag 
maneuver,  curve  shown  the  yaw  angle  and  rudder  angle  in 
degree  as  a  function  of  time,  for  this  study  the  same  pro- 
gram was  used,  but  set  DR  in  dynamic  section: 


DR=  -0.06  98  4 (RAMP (0.0) -RAMP (5.0) ) +0.04 6 56 (RAMP (40.0) . . . 

-RAMP(65.0)-RAMP(14  5.0)+RAMP(170.0)+RAMP(2  50.0)  

-RAMP (265.0) -RAMP (345.0)+  RAMP (425.0) +RAMP (440.0)) 

and  use  prepare  statement  prepare  X,  YAWD  and  prepare  X, 
D00  (YAWD  =  YAW*  57.273,  DOO  =  DR*  57.273). 
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B.   INCLUDED  NONLINEAR  TERMS 

The  computer  program  2  is  the  same  as  the  computer 
program  1,  but  nonlinear  terms  are  added  by  setting: 
NA   =  NA1  +  NA2  +  NA3  +  NA4 


where 


where 


where 


NA1  =  -1(X   q2  +  X   r2  +  X   rp) 
qq      rr      rp  c 

NA2  =  - (mvr  +  X   vr  +  X   wq  +  mwq) 

vr      wq  ^      ^ 

NA3  =  -(X   v2)/l 
vv   " 

NA4  =  -(Alu2  +  A2uu   +  A3u  2)/l 

NB   =  NB1  +  NB2  +  NB3 


NB1  =  -1Y   pq 

pq  M 

NB2  =  -(Y  wp+Y  I  I v I r I +  mur  +  mwp) 
wp  c  v  I  r  |  '  '  tr* 

NB3  =  -(Y   wv  +  Y,  |  |v]v)/l 
wv       j  v  |  v  '  ' 

NF   =  NF1  +  NF2  +  NF3 


NF1  =  -N   pq  +  (Iy-lx)  pq 

NF2  =  -(N   wp  +  N,  ,  |v|r)/l 
wp  ^     j  v  j  r '  '  ' 


NF3  =  -(N   wv  +  N 


wv 


|v|v 


v|v)/l' 


Again  setting  terms  that  include  w,  p,  and  q  (heave,  roll 
and  pitch)  equal  to  zero,  and  set  number  of  known  terms 
from  Table  II  into  section  1  of  CSMP  program  (unknov/n  co- 
efficients set  equal  to  zero)  . 


44 


TABLE  II 


Hydrodynamic  Coefficients  of  Ship  "D"  for  Nonlinear 
Terms  (non-dimensional) 


X    =   0.00005 
rr 

X    =   0.00241 
vr 

X    =  -0.00341 
vv 

Y.  |  ,=  -0.0416 
v|v| 

N  |  i=  -0.01002 
V|  vj 

PROPULSION  RATIO  £  T] 

n  >_  0.45    -l.o  <  n  <  0.45    n  <  i.o 

Al      -0.00004      -0.00032  -0.00117 

A2      -0.00035       0.00070  -0.00100 

A3       0.00099       0  -0.00085 

Fig.  8  -  Fig.  12  are  the  same  results  as  Fig.  3  -  Fig.  7 
The  results  of  computer  program  2  are  more  accurate  than 
the  computer  program  1,  when  compared  with  free  running 
model  test  of  NSRDC  [Ref.  6].   Fig.  13  and  Fig.  14  when 
studying  the  stability  of  the  ship  by  applying  a  force 
moment  to  the  ship  in  computer  program,  set 

NF1   0.0001* (STEP(IO.O) -STEP (10. 01) )  (can  use  NF1  be- 
1  ause  NFl  in  this  program  equal  zero) 

Study  direction  stability  of  the  ship  by  plotting 
direction  of  the  ship  (used  advance  VS.  transfer  ship 
length) and  check  heading  angle  of  the  ship  by  plotting 
YAW  VS.  TIME. 
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IV.   CONCLUSIONS 

The  equation  of  motion  of  surface  ship  and  computer 
program  developed  here  including  all  of  six  degrees  of 
freedom,  but  the  study  in  III  concerns  only  three  degrees 
of  freedom  (surge,  sway,  yaw)  because  hydrodynamic  coeffi- 
cients are  not  available;  when  the  state  of  the  art 
reaches  the  stage  in  which  hydrodynamics  coefficients  are 
available,  this  computer  program  can  be  used  in  all  six 
degrees  of  freedom. 

Some  results  from  III  are  not  too  perfect  because  the 
lack  of  some  constants  and  coefficients  such  as  the  value 
of  mass  (m) ,  initial  velocity  (ADOTO) ,    command  speed  (UC) , 
etc.   But  for  study  can  adjust  from  curve  for  the  model 
test  [Ref .  6] . 

This  computer  program  did  not  include  some  external 
effects  such  as  effects  of  wave  and  wind,  but  these  effects 
could  be  included  in  the  program  by  adding  terms  to  the  IF 
equation . 

The  following  implementations  are  suggested  for  the 
future  work. 

A.  Study  all  six  degrees  of  freedom. 

B.  Study  for  the  effects  of  waves  and  wind. 

C.  Study  for  control  of  the  velocity  and  direction 
of  the  ship  (by  use  of  "MACROS',  "PROCEDURE"  or 
subprogram  in  CSMP) . 
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